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Abstract 

We generalize the following classical result of Fubini [12] for pseudo-Riemannian metrics: 
if three essentially different metrics on M n - 3 share the same unparametrized geodesies, and 
two of them (say, g and g) are strictly nonproportional (i.e., the minimal polynomial of 
g ia g a j coincides with the characteristic polynomial) at least at one point, then they have 
constant curvature. 



1 Introduction 



1.1 Definitions and results 

Two Riemannian or pseudo-Riemannian metrics on the same manifold are said to be geodesically 
equivalent, if they have the same geodesies considered as unparametrized curves. Two metrics g 
and g are strictly nonproportional at a point x G M, if the minimal polynomial of the (l,l)-tensor 
Gj := g tk Qkj coincides with the characteristic polynomial. If one of the metrics is Riemannian, 
strict nonproportionality means that all eigenvalues of G have multiplicity 1. 

The main result of the present paper is the following 

Theorem 1. Let g, g and g be three geodesically equivalent metrics on a connected manifold M n 
of dimension n > 3. Suppose there exists a point at which g and g are strictly nonproportional, 
and a point at which g, g and g are linearly independent. Then, the metrics g,g and g have 
constant curvature. 



If the metrics are Riemannian, the local version of Theorem [T] was proved by Fubini in [T2l 
[13] . The proof of Fubini is short and elegant, but unfortunately, does not work in the pseudo- 
Riemannian case. More precisely, Fubuni's proof is based on the following classical result: 

*MSC: 53A20, 53B21, 53C22, 53C50, 53D25, 70G45, 70H06, 70H33, 58J60 



1 



Theorem 2 (Levi-Civita [20] ). Assume g and g are Riemannian geodesically equivalent metrics 
that are strictly nonproportional at a point p. Then, in a neighborhood of this point there exists 
a coordinate system ( ) such that the metrics g, g are as follows: 



"4 



E 

i=i 



E 



i=l 



JJ^xO-X^x,-)) 



dx? 



dx„- 



(1) 



(2) 



where is a positive function of the variable Xi only. Moreover, every two metrics of this form 
are geodesically equivalent. 



Fubini calculated the riemannian curvature Rij t km of the metric ([I]) and observed that the com- 
ponent Rij^km is non-zero only if = (k,m) or if = (m,k). Nowadays, such curvature 
tensors are called diagonal, see for example [57]. Then he showed that, unless the curvature is 
constant, the curvature tensor determines the coordinates lines (xi, x n ) uniquely Thus, the 
metrics g, g and g must be simultaneously diagonalisable in a certain coordinate system and a 
short analysis shows that they are linearly dependent at every point. 

Our proof, which is also valid in the pseudo-Riemannian setting, is as follows. We study the 
partial differential equation ([5} which is responsible for the fact that the metric g is geodesically 
equivalent to g. The unknown functions in this equation are the components of a (0, 2)-tensor ay 
canonically constructed by metric g. Then, we find the integrability conditions for this equation. 
This is a system of linear equations on a and the Hessian of A := \a a pg a ^; the coefficients in these 
equations are algebraic expressions of entries of the curvature tensor and the metric. We show 
that if the system has two solutions, then for certain K £ R the components of the curvature 
must satisfy the condition afZ ajk i + a"Z aifc/ = 0, where Z ijM = R ijM - K ■ (gug jk - g%k9ji)- 
Then, it is an easy exercise in linear algebra to show that if in addition the metrics g and g 
are strictly nonproportional, then this condition on the curvature implies that the curvature is 
actually constant. Then, the curvature of the metrics g and g is constant as well by the Beltrami 
Theorem. 

Remark 1. We emphasize that the essential part of our proof is, in fact, the analysis of algebraic 
properties of the integrability conditions for the "geodesic equivalence equation" ([5]). Algebraic 
aspects of the integrability conditions will also be clarified in Section [3] where we discuss an unex- 
pected relationship between geodesically equivalent metrics and the so-called sectional operators 
on semisimple Lie algebras. As an application of this observation, we give an alternative, pure 
algebraic proof of the (local version of the) Fubini theorem, and we believe that this idea might 
be useful in a wider context. 
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Remark 2. All assumptions in Theorem [T] are important. Indeed: 

• if the dimension n is 2, counterexamples were constructed in [16], see also [8]; 

• if the metrics are not strictly nonproportional, counterexamples can be found in §9 of [IE], 
see also [Ml HH1 ED]; 

• if we omit the assumption that g, g, and g are linearly independent, then we can take g = g, 
and the Levi-Civita metrics ([I]) and (T5]) for generic X{ give a counterexample to Theorem HJ 

1.2 Motivation 

The first motivation, which was the reason why Fubini studied this question, came from the study 
of projective vector fields of pseudo-Riemannian metrics. Recall that a vector field v is projective 
with respect to the metric g if its (local) flow takes geodesies to geodesies. Projective and affine 
vector fields are very classical objects of investigation: both Lie [21] and Schouten [H] explicitly 
formulated the problem of constructing all metrics admitting one or sufficiently many projective 
vector fields. As a direct corollary of Theorem [1] we obtain 

Corollary 1. Suppose that g and g are geodesically equivalent metrics of a nonconstant curvature 
on a connected manifold M n of dimension n > 3. Suppose there exists a point such that g and 
g are strictly nonproportional. Then, the dimension of the space of projective vector fields minus 
the dimension of the space of homothety vector fields is at most one. 

Recall that the vector field wis a homothety vector field, if L v g = const • g, where L v is the Lie 
derivative. We allow the case const = 0, so Killing vector fields are also homothety vector fields. 

Remark 3. In dimension 2, Corollary [1] is wrong. Indeed, Darboux-superintegrable metrics admit 
three projective vector fields and, as a rule, only one homothety vector field, which is a Killing 
vector field. The definition of Darboux-superintegrable metrics and a description of their pro- 
jective vector fields can be found in [8]. Moreover, it is possible to show that there are no other 
counterexamples to the 2-dimensional version of Corollary [TJ This is a very nontrivial statement 
which follows from the results of [8j [38], where all 2-dimensional metrics admitting projective 
vector fields were constructed. 

Moreover, it is possible to use Theorem [1] in order to describe all projective vector fields of a 
metric g under the assumptions that there exists a metric g which is geodesically equivalent to 
g and strictly nonproportional to g. Indeed, in this case the system of PDE on v is a system of 
ODE's, which is much easier to analyse than a system of PDE's, seefllElEliniEIlElESlEZlEH], 
where this ODE-system was obtained and completely solved in particular cases. 

One more motivation comes from the theory of superintegrable systems. Recall that a metric is 
called superintegrable, if the number of independent integrals of special form is greater than the 
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dimension of the manifold. There are different possibilities for the special form of integrals; de- 
facto the most standard special form of the integrals is the so-called Benenti-integrals, which are 
essentially the same as geodesically equivalent strictly nonproportional metrics, see [4J, [TJ, [TTj . The 
results of the present paper show that Benenti-superintegrable metrics of nonconstant curvature 
cannot exist, which was a folkloric conjecture. 

1.3 History 

The theory of projective transformations has a long and fascinating history. The first non- 
trivial examples of projective transformations were discovered by Lagrange [T9l [35] . Geodesically 
equivalent metrics were studied by Beltrami, Darboux, Levi-Civita, Painleve, Eisenhart, Weyl 
and other classics. One can find more historical details in the surveys [31 HDJ [9] and in the 
introduction to the papers [371 EQl [291 ESI [Ml EHl [28l [27]. 

The pseudo-Riemannian version of Fubini's Theorem was investigated by the Kazan school of 
geometry, in particular by A. Aminova and her collaborators. In particular, they tried to mimic 
the Fubini proof for pseudo-Riemannian metrics. They used the description of geodesically 
equivalent pseudo-Riemannian metrics obtained by Aminova [2J, which is a generalization of 
Levi-Civita's Theorem [2j Unfortunately, the description by Aminova is very complicated, which 
makes this programm to be computationally very hard. Besides, there are infinitely many different 
types of normal forms for pseudo-Riemannian metrics (depending on the Jordan form of the 
tensor G), and each of these types require a separate inverstigation. According to [HE], they 
proved Theorem [T] in dimensions up to 6, though we did not find the place where the proof is 
written. 
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2 Proof 



2.1 Schema of the proof 

We first show (Section I2.3j) that if geodesically equivalent metrics are strictly non-proportional 
at least at one point, than it is so at almost every point. Then, we will show (Section 12.41) that if 
g, g, and g are linearly dependent at every point of a neighborhood, they are linearly dependent 
at every point of the manifold. 

Combining these two observation, we see that if the assumptions of Theorem [I] are fulfilled at one 
point, they are fulfilled at almost every point, i.e., though the theorem is global, it is sufficient 
to prove it locally. This will be done in Sections 12.51 12.61 



2.2 What we will use in the proof 

There are many tensor reformulations of the condition "metrics g and g are geodesically equiva- 
lent". In our paper we will use the following one, which was suggested by Sinjukov [37], see also 
[3 HD]: given two metrics g and g, consider the (0,2)— tensor 



det(fl) 



det(g) 



i 

,a(3 



gi a g a p gjp, (3) 



and the function 

A := laaptf* (4) 

where and g alS and g a/3 denote the tensors dual to <7y and respectively, i.e., g ia g a j = 5j and 
9 ia 9aj = S). 

Theorem 3. JJ1\ [7]/ The metrics g and g are geodesically equivalent, if and only if 

0>ij,k = ^,igjk + \jgik, (5) 

where the covariant derivative is taken with respect to the (Levi-Civita connection of the) metric 
9- 

We will also use the following connection between geodesically equivalent metrics and integrable 
geodesic flows due to 



Let Oj be as in (j3J) , with one index lifted by g (so now a is a (1, 1)— tensor, self-adjoint with 
respect to g). Consider the family S t of (1, l)-tensors 

S t = det(a - t Id) (a - t Id)" 1 , tel. (6) 
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Remark 4. Although (a — t Id) - is not defined for t lying in the spectrum of a, the tensor St is 
well-defined for every t. Actually, the matrix of St is the comatrix of the matrix of a — t Id. In 
particular, St is a polynomial in t of degree n — 1 with coefficients being (l,l)-tensors. 

We will identify the tangent and cotangent bundles of M n by g. This identification allows us to 
transfer the natural Poisson structure from T*M n to TM n . 

Theorem 4 (pU [52J, [53]). If a is a solution of ([5j) ; i/ien /or even/ feR t/ie function 

I t :TM n ^R, I t (0 = 9(S t (0,0 (7) 
is an integral for the geodesic flow of g. 

Recall that a function is an integral of the geodesic flow, if it is constant along the trajectories, 
i.e., along the curves on the tangent bundle of the form (^(t), 7(t)), where 7 is a geodesic and 
j(t) is its velocity vector at a point t. 

We will also use the following statement, whose Riemannian version was obtained in [24J, and 
pseudo-Riemannian generalisation is due to Topalov |51j. As in the introduction, we denote by 
G the (1, 1) -tensor G) := g ik (jkj- 

Theorem 5 (follows from Theorem 2 and Section III B of [51]). Suppose the degree of the 
minimal polynomial of G is r at every point of the neighbourhood U(p). Then, for arbitrary 
mutually different t\,...,t r+ i G R the functions I tl ,...,It r are functionally independent, and the 
function It r+1 is a linear combination of the functions lt x ,---,lt r (the coefficients of the linear 
combination are constant). 

Recall that functions fx,...,f r are functionally independent, if their differentials are linearly in- 
dependent almost everywhere. 

In other words, the number of independent integrals among I t is the degree of the minimal poly- 
nomial. In particular, if the metrics are strictly nonproportional at a point p (which immediately 
implies that they are strictly non-proportional at every point in a small neighbourhood U(p)), 
then the differentials of the integrals are linearly independent at almost every point of TU(p). 

2.3 If geodesically equivalent metrics are strictly nonproportional at 
one point, then they are strictly nonproportional at almost every 
point. 

Let geodesically equivalent metrics g and g be strictly nonproportional at p. Consider a geodesic 
7 passing through p := 7(0). Let us show that every point q := j(r) of the geodesic has a 
neighborhood U (q) such that at almost every point of the neighborhood the metrics are strictly 
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nonproportional. Since every point can be reached by a finite sequence of geodesies, and since 
the condition "the minimal polynomial of G has degree n" is an open condition, this will prove 
the statement formulated in the title of this section. 

As we recalled in Section 12.21 for almost every point of TU (p) the differentials of the integrals 
I tl ,...,It„ are linearly independent. Take a sequence of points (pk, v k) £ TU(p) converging to 
(7(0), 7(0)) such that dl tl (pk, Vk), ■ ■ ■ , d!t n {Pki v k) ar e linearly independent. Consider the sequence 
of the geodesies 7^ such that 7^(0) = pk and 7&(0) = v^. For sufficiently large k, the geodesies 
exist up to time r and the sequence ( / jk( T ), 7fc( T )) converges to (7(1"), 7 (r)). 

Since the integrals are preserved by the geodesic flow, the differentials of the integrals are also 
preserved by the geodesic flow. Then, at the points (7fc(r), jk(r)) the differentials of I tl , . 
are linearly independent. By Theorem 0, in an arbitrary small neighborhood of 7&(t) there exists 
a point qk such that the metrics are strictly nonproportional at implying the claim, □ 

2.4 If gik, gkj, and g^j are linearly dependent at every point of a 
neighborhood U, then they are linearly dependent at every point 
of M. 

Within this section we assume that g, g, g are geodesically equivalent metrics on a connected 
manifold M of dimension n > 3. We consider the tensor and the function A given by fl3El) . 
The same objects for the pair of metrics g, g will be denoted by the capital letters A and A, i.e., 

1 • g ia g^g jP , A := \A f #<f# (8) 
We will first prove the following statement (essentially due to Weyl |56j). 

Lemma 1. Suppose a and A are solutions of (JS]). Assume a = C ■ A, where C is a function. 
Then, C is a constant. 

Proof. Our proof is different from the proof of Weyl and is based on the ideas developed in 
[53]. Note that in the proof we use only the fact that the dimension is greater than one, i.e., it 
works in dimension 2 as well. 

Consider the integrals Jo from (J7J) constructed by a (we keep the notation Iq for it) and by g 
(we denote it by Xq). If a = C ■ A, then the integrals Iq and Zo are proportional as well, direct 
calculations show that /o(£) — ±C n_1 -X (^). Since the functions Iq and X are constant along 
every trajectory of the geodesic flow, the coefficient of proportionality of these functions is also 
constant along every trajectory of the geodesic flow implying that it is constant everywhere. □ 

Now let us assume that at every point of the neighborhood U the tensors g, g, and g are linearly 
dependent. Then, for certain functions c, d the tensors g,a,A satisfy (probably in a smaller 



A 



det(g) 



det(g) 



7 



neighbourhood U' C U; without loss of generality we can think that U' coincides with U.) 

(Ijj Cflij ■ d.\, r (9) 

We will show that the functions c, d are actually constants. 

Differentiating ([9]) and substituting ([5]) and its analog for the solution A, we obtain 

^,i9jk + X,j9ik = c k gij + dA ti g jk + dA d g ik + d >k Ay, (10) 
which is evidently equivalent to 

ngjk + Tjg ik = c^gij + d^A^, (ll) 

where Tj = — dA^. We see that for every fixed k the left-hand side of ( fTTT) is a symmetric 
matrix of the form TiVj + TjVi. If c& is not proportional to d jk , this will imply that g^ also is of the 
form T^j + TjVi, which contradicts the non-degeneracy of g. If c jk = / • d tk , then the coefficient / 
of the proportionality should be a constant implying d = const • c and = c (g^ + const ■ A^). 
Since the equation (jSJ) is linear and since g^ and Ay are its solutions, their sum g^ + const • Ay is 
also solution of (JSJ). Then, Lemma [T] implies that c is constant. Thus, c and d are constant in a 
neighbourhood U. Since the equation (jSJ) is linear and of finite type, see [TO], linear dependence 
of solutions in a neighbourhood implies linear dependence of the solutions everywhere. □ 

Remark 5. Though we used that the dimension of the manifold is at least three, the statement 
is true in dimension two as well provided the curvature of g is not constant, see [15] . 



2.5 Main step of the proof of Theorem [T] 

As we explained in Section 12.11 in view of Sections 12. 3[ \2A\ it is sufficient to prove Theorem 
locally only. 

Assume the metrics g, g, and g are linearly independent at a point p. Assume the metrics g, g are 
strictly nonproportional at p. Then, for a certain neighborhood U(p) the metrics g,g, and g are 
linearly independent, and the metrics g,g are strictly nonproportional at every point of U(p). 

Then, the solution a corresponding to g, the solution A corresponding to g, and the metric g are 
linearly independent at every point of U(p), and the minimal polynomial of the (l,l)-tensor a* 
has degree n. 

As in Section 12.31 we denote by A the function (jlj) constructed by A. 

Let us first consider the case when A is constant. In this case, the equation (jSJ) implies that ay is 
covariant constant. Since a 1 - is self-adjoint with respect to g, at every point there exists a basis 
(bi(p), ...,b n (p)) such that the matrix of g is diagonal with ±1 on the diagonal and the matrix of 
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cij is Jordan-matrix. The basis is unique up to signs of the vectors, so locally b\ can be taken to 
be smooth vector fields on the manifold. Since the vectors of the basis are invariantly constructed 
by two covariantly-constant object, they are covariantly constant as well. Then, the metric g is 
flat implying that g and g have constant curvature. 

In what follows we assume that A is not constant. 

Theorem CD is an easy consequence of the following two lemmas. 

Lemma 2. Assume that a and A are solutions of (jSJ) such that a,A,g are linearly independent. 
Then, at every point there exists K such that the tensor Z^i '■= Rij,ki~ K ' ' (du9jk — 9ik9ji) satisfy 
the condition 

0, (12) 

We see that by construction the tensor Z^i is skew-symmetric with respect to the first two 
indexes. Note that in Lemma [2] we did not assume that the minimal polynomial of a has degree 
n, i.e., the lemma is valid in a slightly more general setting than we need it. 

Remark 6. Tensor Z plays important role in the theory of geodesically equvialent metrics; it 
appears quite naturally in the investigation of special metrics such that Einstein metrics, pseudo- 
symmetric metrics, Kahler metrics, see for example [bH [15] for details. 

The next lemma shows that if in addition the minimal polynomial of a has degree n, then the 
condition (|2o^) implies Z^i = 0. Then, the curvature of g is constant. By Beltrami's Theorem 
(see for example [35]) or the original papers [3] and [IS]), the metrics g and g have constant 
curvature as well which proves Theorem [TJ 

Moreover, we will not use the indexes k and / in the proof, so Lemma [3J is the matrix reformulation 
of the condition afZ a j + a^Z ai for arbitrary (1, 1)— tensor a and a skew-symmetric (0,2)-tensor 
Z. 

Lemma 3. Let Z , a be n x n-matrices such that Z is skew- symmetric. Assume the minimal 
polynomial of a has degree n. If Za — o^Z = 0, then Z = 0. 

The proof of Lemma [3] is an easy exercise in linear algebra and will be left to the reader (it is 
a kind of problem such that it is easier to prove than to understand the proof.) We recommend 
to consider the coordinates such that the matrix a is in Jordan form, and then to calculate the 
matrix Za. One immediately sees that if the matrix Za is symmetric, then Z = 0. □ 

2.6 Proof of Lemma [2]. 

The proof is straightforward tensor calculations. The geometry behind the calculation could 
be understood with the help of [10]. There, the equations © were written in the projectively 
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invariant form (so that the equations for g and for the geodesically equivalent metric g are the 
same). The prolongation of the equations was written as a connection on the projective tractor 
bundle, and the curvature of the connection was calculated. Its first part related to (the analog of 
the tensor) a is a trace-free object. Moreover, for most objects, the trace-free part of the covariant 
derivative coincides with the trace-free part of the corresponding derivative on the tractor bundle. 

We will see that in the proof we consider the integrability conditions for the equations (151) , and 
then "artificially" write all objects in trace-free form. At the end we obtain the required equation 

(USD- 

Note that in the proof we will essentially use the symmetries of the Riemannian curvature tensor, 
which have sense only if the afline connection is the Levi-Civita connection of a metric. That 
means, it is important for us that among the solutions of the projective-invariant analog of a 
there is a non-degenerate solution. 

In Section [3], we give another proof of Lemma [2] which is, in fact, pure algebraic and based on 
some ideas from the theory of integrable Hamiltonian systems on Lie algebras. 

Proof of Lemma [21 Integrability conditions for the equation ([5]) are (we use the standard fact 
that dij^i — dijik = a ia R^ kl + a a jR^ kl for any (0, 2)— tensor a^) 



a iaR'jkl + a ajR C ikl — \li9jk + ^,lj9ik ~ \ki9jl ~ ^,kj9il- (13) 

The same is true for the other solution A 

AiaR'jkl + AajRfkl = \li9jk + ^,lj9ik ~ ^,ki9jl ~~ ^,kj9il- (14) 

Starting from this point, the proof is purely algebraic: the statement of Lemma [2] is an algebraic 
corollary of ffTBl and ffl~4l) . We emphasize this once again in the next section by giving another 
version of the proof in the Lie-algebraic language. 

Now let us multiply the equation (1131) by A\ and sum over I. After renaming indexes, we obtain 
a iaR'jk^Af + a a jRf k/3 Ai = X^aiA^gjk + X^A^g ik — X^Aji — X^jAu. (15) 

Using the symmetry of the Riemann tensor we obtain a^R^^ — a iRaj,kf3^i = afRpkjaAi = 
a^ApiR^-. Substituting this in (TToT) . we get 



a ?ApiR kia + a"A/3iR k j a — X^iA^gjk + X^A^g,-^ — X^iAji — X^jAu. (16) 

Let us now symmetrise (1161) by I, k to obtain 

af (A P iRl ja + ApkRf^ + a a j (A pk Rf ia + A P iRP ki ^j ^ 
= X^aiAfgjk + X >a jAfgik — X^i.Aji — X^jAu + X ya iA k gji + X )Ct jA k gu — X t uAjk — XjjAif.- 
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We see that the components staying in brackets are the left-hand side of the equation (jT3j) with 
other indexes. Substituting ()14p in the first term of the left-hand side we obtain 



O-i {^-piRkja + ApkRfjc^j = af (-^,al9jk + ^-,ak9jl ~ — Ajkgia) 

— a i^,al9jk + a i-h-,ak9jl ~ A-,jl a ik — \jk a il- 

Similarly, the second term of the left-hand side is 

a" (As^L + A pkRuo) = a] (^,ai9ik + A ,ak9il — \iWka ~~ \ik9la) 

= a^A ja ig uk + a'jA^kgu — A^a^ — A jik dji. 

Substituting (I18II19P in (ITT]) , we obtain 

a>iA t ai9jk + afA^gji — Ajidik — Aj k du + d^A )a ig ik + d^A^ ak gu — A t uaj k — A )ik dji 

= X^A^g^ + ^, a jAfgik — X^iAji — X )k jAu + X^A^gji + X^A^gu — X t uAj k — X t ijA ik . 

Collecting the terms by g, we see that (1201 can be written as 



(19) 



(20) 



(a?A >al - X, ai Af) g jk + (a?A >ak - X, ai Al) gjl + (ajA >al - X, aj Af) g ik + (afA >ak - X, aj A^) g u 
= Aj^ajfc + Ajkdu + A^icijk + A^kdji - X^iAji — X^jAn — X.uAjk — XjjAik. 

(21) 

We denote Tu := afA^i — AfX >a i- In this notation, the equation (12~T1) is 

TiWjk + Tikgjl + TjWik + Tjkgil ^ ^2) 

= Ajidik + Ajkdu + A^djk + A. ik dji — X t kiAji — X t kjAu — X t uAj k — XjjA ik . 

Let us show that r is symmetric. Multiplying with gi h and contracting with respect to j, k, we 
obtain 



(n + 2)r a + (r jk g jk ) g u = A >al a? + (g jk A, jk ) d U + (d jk gi k ) + A :ia df 

-X, ai Af - {X, kj g k i) An - X M {A jk g^ k ) - X M A? 



(23) 



We see that the right-hand side is symmetric w.r.t. i, I. Then, so should be the left-hand-side 
implying r»j = Tu. We also see that the sum of the first, fourth, fifth and last terms of the 
right-hand side is Tu + tu = 2r^. Then, the equation f!23p is equivalent to 

ra = -(- {r jk g jk ) gu + {g jk Aj k ) a a + A )U (d jk g jk ) - (X, kj g kj ) A u - X M {A jk g jk )) . (24) 
Now we return to the equation ( 1221) . We alternate the equation with respect to j, k. 

T ikgjl + T jWik — T ijgkl — T kWij ^5) 

= Ajidi k + A )ik dji — X, k iAji — XjjA ik — A )k idij — A^jd k i + XjiA k i + Xj k Aij. 
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Let us now rename % <-> k in (|25|) . We obtain 



Adding fl22|) with (126]) and dividing by 2 for cosmetic reasons, we obtain 



Tikdjl + 7"j7<?ifc — Ajid ik + A jifc a — A^jA^ — \jjA ik . 

Substituting the expression for r from ( l24l) . we obtain 



(W") 9ik + (g al3 A a p) d ik + K lk (d a pg a P) - (X^g^) A lk - \ ki (A a pg^\) g jt 
+ ~ (- (r a ^ Q/3 ) g jt + {g aP A,*p) a jt + A s (aapgrt) - (\ a pg aP ) A fl - A, y (A,^)) ^ 
= Ajidi k + A^kdji — X,uAji — XjjAik- 

Denoting (r^g 130 ) by r, (A,^^) by /i, (A jj3a ^ a ) by M, and using that (a a pg a P) = 2 A, and 
(A Qi g<7°^) = 2A, one obtains 

\ {-TQik + Md ik + 2A iik X - fiA ik - 2A ifci A) g jt 
+ \ (-rgji + Maji + 2Aj,A - nA iX - 2A, y A) g ik 
= Ajid ik + A )ik d ji — X tk iAji — X t ijA ik , 

Combining the terms, we obtain 



i-lg* + 2A, ik X) f - (\g lk + X M A) f + (-Zg fl - 2X tlj A) * + (-Zg fl + 2A Jl X) 
= {--Mgji + Aji) d ik + [A jik - ^Mg ik ) a jt - (X M - ^fig*) A jt - (X :tj - ^gji) A ik . 

(27) 

Now, let us calculate r: we multiply ffMb by g d and sum over i, I. After dividing by 2, we obtain 

r = - (MA - wA) . 

n 

Substituting in fl27l) . we obtain 



[—^9ik + A,ifc) ^7 - ("nto + ^,fc») ^ndjl + (n9jl ~ ^9ik + \—^9jl + ^9ik 

{-\Mgji + Aji) d ik + (A tik - \Mg ik ) aji - (X tki - \ngik) A jt - (X^ - ^figji) A ik , 
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which is equivalent to 



= (A,« - M (Afl - f 9jl ) + (A, M - (A, - ^) ^ 

Denoting S ife := (A ftt - ^^«) , &tt := - d/i := (o^j - ^foj), D jt := (A ik - ^g ik ), 

we see that (I28I) is equivalent to 

Bikdji + Bjidik = bikDji + bjiDik- (29) 

Now, it is easy to see that the condition B a dp + Bpd a = b a Dp + bpD a for dimensions > 3 implies 
that B is proportional to b and D to d, or that B is proportional to D and 6 to d. We see that the 
condition (|29p is essentially the same as this condition, the role of a and f3 play the multi-indexes 
ik and jl. Since g, a, and A are linearly independent, D can not be proportional to d. Thus, b 
is proportional to d which imply that X^j is a linear combination of Oy and g^. The coefficients 
of the linear combination are not important for us, let 

\ij = P- 9ij + K ■ Ojj 
Substituting this equation in ffTB"]) . we obtain 



a iaZf kl + a a jZ^ kl — 0, (30) 
Where Zj w = R % - kl — K ■ (8\gj k — b k 9ji)- Clearly, the equation (1311 is equivalent to the equation 

We see that by construction the tensor Z^i is skew-symmetric with respect to the first two 
indexes (actually, the tensor Z has the same symmetries as the curvature tensor, i.e., for example, 
skew-symmetric with respect to the last two indexes as well, but we will need only the first two 
indexes). Thus, Lemma [2] and, therefore, Theorem [1] are proved, □ 



3 Fubini theorem and sectional operators 
on semisimple Lie algebras 

In this section we discuss an unexpected and remarkable relationship between geodesically equiv- 
alent metrics and some special operators on semisimple Lie algebras which appeared in the theory 
of integrable systems. 

We start with a brief overview on (one special type of) integrable Euler equations on semisimple 
Lie algebras (see QU ED H21 03] for details). 
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Let g be a semisimple Lie algebra, R : g — > g an operator symmetric with respect to the Killing 
form ( , ) on g. The differential equation 

x—[R(x),x], x G g, (31) 

is Hamiltonian on g with respect to the standard Lie-Poisson structure and called the Euler 
equation related to the Hamiltonian function H(x) = ^(R(x),x). 

A classical, interesting and extremely difficult problem is to find those operators R : g — > g for 
which the system fl3Tl) is completely integrable. 

One of such operators was discovered by S. Manakov in [22] and his idea then led to an elegant 
general construction developed by A. Mischenko and A. Fomenko [JT] and called argument shift 
method. This construction in brief can be presented as follows. 

Assume that R : g — > g satisfies the following identity 

[R(x),a] = [x,b], x G g, (32) 

for certain a, b G g, a ^ 0. Then the following statement holds 

Theorem 6. [JT] Let R : g — > g be symmetric and satisfy ( 1321) . Then 

1 ) the system (l3Tj) admits the following Lax representation with a parameter: 

— (x + Aa) = [R(x) + Xb, x + Aa]; 

(Jbb 

2) the functions f(x + Aa), where f : g — > R zs an invariant of the adjoint representation, are 
first integrals of (J3"I|) /or any A G R and, moreover, these integrals commute; 

3) if a E g is regular, then (13T1) zs completely integrable. 

This construction has a very important particular case. If the Lie algebra g admits a ^-grading, 
i.e., a decomposition g = (} + o (direct sum of subspaces) such that [fj, fj] C f), [f), t)] C D, [o, o] C f), 
then we may consider R : f) — > ^ satisfying (1521 with a, 6 G 0, and Theorem still holds if we 
replace g by [). 

The most important example for applications (in particular, in the theory of integrable tops) is 
g = sZ(n,R), t) = so(n, R), with a and 6 symmetric matrices. This is the situation that was 
studied in the pioneering work by S. Manakov |22j leading to integrability of the Euler equations 
of n-dimensional rigid body dynamics. 

From the algebraic point of view, the above construction still makes sense if we replace so(n) by 
so(p, q) and assume a, 6 to be symmetric operators with respect to the corresponding indefinite 
form g. Moreover, if we complexify our considerations we do not even notice any difference. 
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However, to indicate the presence (but not influence) of the bilinear form g, we shall denote the 
space of g-symmetric operators by Sym(g), and the Lie algebra of (^-skew-symmetric operators 
by so(g). 

Definition 1. We shall say that R : so(g) — > so(g) is a sectional operator associated with a,b G 
Sym(g), if R is symmetric with respect to the Killing form and the following identity holds: 

[R(x),a] = [x,b], for all x G so(g). (33) 

We follow the terminology introduced by Fomenko and Trofimov in [5U [11] where they studied 
various generalizations of such operators. Strictly speaking, the above definition is just a partic- 
ular case of a more general construction. The term "sectional" was motivated by the following 
reason. The identities ( 1321) and ( 133]) suggest that one may represent R as ad~ 1 ad; ) , but in general 
we cannot do so because ad a , as a rule, is not invertible. That is why the operator R splits into 
parts each of which acts independently on its own subspace (section). 

A surprising relationship between sectional operators and geodesically equivalent metrics is ex- 
plained by the following observation. Notice, first of all, that due to its algebraic symmetries 
(skew-symmetry with respect to i,j and k, I and symmetry with respect to permutation of pairs 
(ij) and (hi)), the Riemann curvature tensor Rij^i can be naturally considered as a symmetric 
operator R : so(g) — > so(g) (strictly speaking we need to raise indices i and k by means of g to 
get the tensor of the form R 1 -^). Having this interpretation of R in mind, we immediately obtain 

Theorem 7. Let g and g be geodesically equivalent nonproportional metrics, then the Riemann 
curvature tensor R\ of the metric g is a sectional operator in the sense of Definition CD More 
precisely, 

[R(x), a] = [x, b] 

where a is the g-symmetric operator associated with the form aij defined by (3), and b is the 
g-symmetric operator associated with the form 2X ^ (Hessian oftva). 

The proof of this statement is just the observation that fl33|) is a translation of the compatibility 
condition ([TBI into Lie-algebraic language. 

Before discussing the proof of the Fubini theorem in this "new" language, we make some remarks 
which could also be useful. 

Notice, first of all, that in our new notation the condition "curvature is constant" (at a point) 
simply means that R : so(g) — > so(g) is a scalar operator, i.e. R(x) = K ■ x. 

Furthermore, it is a very simple Lie-algebraic fact that fl33|) implies that a and b commute. 
Indeed, ([b,a],x) = (a, [x,b]) = (a, [R(x),a]) = ([a, a], R(x)) = for any x G so(g), so [a, b] = 0. 
In the theory of projectively equivalent metrics this means that the operator a commutes with 
the Hessian of its trace 2A = tr a. This fact is, of course, well known, but the above proof seems 
to be the simplest one. Moreover, if instead of a we substitute any element £ from its centralizer 
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C(a), we obviously get the same conclusion [6, £] = 0, i.e., b lies in the center of the centralizer 
of a. This means, in fact, that 6 is a polynomial of a. 

Finally, if a is regular in the Lie-algebraic sense, i.e. its minimal polynomial coincides with 
the characteristic one, then the operator ad a : so(g) — > Sym(g) has trivial kernel so that the 
sectional operator R (i.e., the curvature tensor!) can be reconstructed from a and b. Namely, 
R(x) = ad~ 1 adf,(x), a well-known formula in the theory of integrable systems on Lie algebras. If 
we take into account the fact that b = P(a) = A n _ia n_1 + A n _ 2 a™~ 2 H — • + A x a + A (polynomial 
of a), then this formula can be rewritten as 

R( x ) = lp(a + tx)\ t=Q . (34) 

Indeed, [P(a + tx) , a + tx] =0 implies 

= j t [P(a + tx),a + tx}\ t=0 = [j t P(a + tx)\ t=0 ,a] + [P(a),x], 

i.e., [jgP(a + £r)| t= o, a] = [x,b]. Since a is regular, we have fl34l . 

This shows, in particular, that the algebraic structure of the curvature tensor can be understood 
in terms of the operator a only. 

Using this language we now give another proof of the tensor part of the Fubini theorem (Lemmas 
2 and 3). 

Assume that we have three geodesically equivalent metrics g, g, and g. Then the Riemann 
curvature tensor R of the metric g satisfies at the same time two identities : 

[R(x),a] = [x,b] and [R(x),A] = [x,B], (35) 

where a) = g ki a lv A) = g kl A iV b) = 2g H \, v B) = 2g ki K M (cf. (TH and (HI). 

From now on, we may forget about the geometrical meaning of a, b, A, B and start thinking 
of them as just certain ^-symmetric operators. In addition, without loss of generality we may 
assume all these operators to be trace free (as, of course, it should be in the semisimple Lie algebra 
sl(n,M) which stands behind this construction). Moreover, we are allowed to complexify all the 
objects so that instead of so(g) and Sym(g) we may simply consider the spaces of symmetric and 
skew-symmetric complex matrices. 

The reformulation of the (algebraic part of) Fubini theorem are the following analogs of Lemmas 
2 and 3 respectively. 

Lemma 4. Let R : so(g) — > so(g) be symmetric and satisfy fl35|) . If a and A are not proportional, 
then b is proportional to a and, therefore, [R(x) — K ■ x,a] =0 for some K £ R. 

Lemma 5. If a is regular, i.e., its minimal polynomial coincides with the characteristic one, then 
the identity [R(x) — K ■ x,a] =0 implies R = K ■ id (i.e., the curvature is constant). 
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Proof of Lemma HI Let y and z be arbitrary ^-symmetric matrices, then [A, y], [a,z] G so(g) 
and we have: 

[R([A,y]),a] = [[A,y],b], [R([a,z]),A] = [[a,z],B]. 
Since R is symmetric with respect to the Killing form ( , ) we have 

{[[A,y],b],z) = {[R([A,y]),a],z) = (R{[A,y)),[a, z)) = ([A,y),R{[a,z])) = 

(y,[R([a,z}),A\) = (y, [[a,z],B]} = ([[B,y],a],z) 
Since z is an arbitrary symmetric matrix, we conclude that 

[[A,y],b] = [[B,y],a]. (36) 

This relation is an analog of ( 1201) . Similarly, [[a,y],5] = [[&, j/],v4]. Using the Jacobi identity, it 
is not hard to see that 

[b, A] = [a, B] 

Rewriting (136]) as 

y(Ba - Ab) + (aB - bA)y = Bya + ayB - by A - Ayb 
and noticing that [b, A] = [a, B] implies Ba — Ab = aB — bA, we get 

yT + Ty = Bya + ayB — by A — Ayb 
where T denotes aB — bA (this is an analog of r from (122]) ). 

This formula can be considered as a relation between two linear operators acting on the space of 
symmetric matrices (the argument of both operators is y e Sym(g)). To get some consequences 
from this identity, we take a kind of its trace. Recall that we consider A, a, B, b, y, T as usual 
symmetric (complex) matrices. 

Instead of y we substitute the symmetric matrix of the form eiV T + vej , where ej and v are 
vector-columns (ei, . . . ,e n is the standard (orthonormal) basis), then apply the result to ej and 
take the sum over i. Here is the result: 

(eiV T + veJ)Tei + T(e i v T + vej)ei = B(eiV T + vej)aei + ... 

ei(Tv, + v(Te i} ej) + Tei(v, + Tvfe, ej) = Be^av, ej) + Bv(aei, ej) + ... 

Using obvious facts from Linear Algebra such as 

y^(Tej,ej) =trT, J^(e i ,e i )=n, e^v, = v , 
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we get 

Tv + tr T ■ v + Tv + n ■ Tv = Bav + tr a ■ Bv + ... 
Taking into account that a,A,b,B are all trace free we have 

((n + 2)T + trT • ld)v = (Ba + aB - Ab - bA)v 

Since u is arbitrary and T = Ba — Ab = aB — bA, we finally get 

nT + tr T ■ Id = 

But this simply means that T = 0. Hence we come to the identity of the form 

Bya + ayB = by A + Ayb. (37) 

It remains to use the following simple statement: if a, b, A, B are symmetric, a^O and fl37|) holds 
for any symmetric y, then either b = K ■ a, or A = K ■ a for some constant Kel. 

By our assumption, a and A are not proportional, so we conclude that b = K ■ a and therefore 
the identity [R(x), a] = [x, b] becomes [R(x) — K ■ x, a] = 0, as needed. 

Notice that ( 1371) and the rest of the proof almost literally repeat ( l29i) and the end of the proof of 
Lemma 2, □ 

Proof of Lemma [5J Let [R(x) — K ■ x, a] = and a be regular. It is a well known algebraic 
fact that the centralizer of a regular matrix a is generated by the powers of a. In particular, 
the centralizer of a consists of (/-symmetric matrices. On the other hand, R(x) — K ■ x is skew- 
symmetric. Thus, R{x) — K ■ x has to be zero for any x, i.e., R = K ■ id, as was to be proved, □ 
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